Momentum Resolved Optical Lattice Modulation Spectroscopy for Bose Hubbard Model 
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We propose a new method of optical lattice modulation spectroscopy for studying the spectral function of 
ultracold bosons in an optical lattice. We show that different features of the single particle spectral function in 
different quantum phases can be obtained by measuring the change in momentum distribution after the modu- 
lation. In the Mott phase, this gives information about the momentun dependent gap to particle-hole excitations 
as well as their spectral weight. In the superfluid phase, one can obtain the spectrum of the gapless Bogoli- 
ubov quasiparticles as well as the gapped amplitude fluctuations. The distinct evolution of the response with 
modulation frequency in the two phases can be used to identify these phases and the quantum phase transition 
separating them. 

PACS numbers: 



Ultracold atoms on optical lattices have opened up new pos- 
sibilities of studying interacting quantum systems with tun- 
able parameters in a controlled environment [ 1 1 . A seminal 
achievement in this field has been the implementation of the 
Bose Hubbard model [2], which is used as the paradigm for 
superfluid-insulator (SI) transitions, and has been studied ex- 
tensively both analytically G]-[6), and numerically f7]|8)- The 
observation of the SI transition |9) in this system is a suc- 
cessful milestone in the study of ultracold atomic systems on 
optical lattices. More complex but related systems like Bose- 
Fermi mixtures [10], spinor Bose condensates ifTTI . and Bose 
Hubbard model with an artificial vector potential created by 
Raman beams lfl2l have also been implemented leading to a 
host of novel phenomena |[T3l4l6l . 

The SI transition in these systems is usually observed 
through time of flight images of boson density which yields 
information about the boson momentum distribution, n^, in 
the equilibrium state [9|. The superfluid phase is then iden- 
tified by sharp peaks of nu at zero momentum (and at corre- 
sponding Bragg vectors), while the Mott phase shows a dif- 
fuse density pattern. The contrast of the peaks (or visibility) 
has traditionally been used to track such SI transitions. This 
is, however, unsatisfactory for a clean distinction between the 
two quantum phases. The main problem with this method, 
as pointed out in earlier works IfTTI . is that has a precur- 
sor peak at k = near the SI transition 0. Thus look- 
ing for a signature of superfluidity involves looking for an 
extra peaky feature on top of the generalized cluster around 
k = which is experimentally challenging. The other ap- 
proach that has been suggested is to measure the local density 
in situ with a very fine resolution and look for a plateau of 
constant density corresponding to the Mott phase lfl8l [191 . 
An alternate direct way of distinguishing between the super- 
fluid and Mott phases is to measure the excitation spectra in 
these two phases. The superfluid phase has gapless spectrum 
with a linear low energy dispersion, while the spectrum in the 
Mott phase is gapped. Indeed the energy transfer due to opti- 
cal lattice modulation ll20l422l has been previously studied as 
a possible probe for the Bose Hubbard model Il23ll24ll . How- 



ever, such a method does not yield momentum-resolved infor- 
mation about the spectral function of the bosons. 

In this letter, we propose a novel method based on optical 
lattice modulation [20-22 1 to obtain the single particle spec- 
tral properties of the Bose-Hubbard model in both the Mott 
and the superfluid phases. (In a broad qualitative sense, our 
proposed technique, as we establish in this work, is some- 
what akin to being the bosonic counterpart of the ARPES 
technique used extensively in solid state electronic materi- 
als in the sense that our method gives direct momentum- and 
frequency-resolved spectral information about the low-lying 
quantum excitations of the interacting system.) Our proposed 
experiment is the following: we turn on an additional weak 
modulating optical lattice at frequency cj for some time be- 
fore turning off all lattice and trap potentials simultaneously 
and measuring the density of the particles after a time of 
flight. For a sufficiently long time of flight, the density in 
real space can be mapped back to the momentum distribu- 
tion of the system after the lattice modulation. We propose 
to compare this momentum distribution to that in the unper- 
turbed case (without the modulation). We note at the outset 
that such a comparison of momentum distributions involves 
looking at the excitations created by the optical lattice modu- 
lation and not measuring equilibrium of the bosons. Hence 
the proposed experiment does not suffer from the problem of 
measuring the equilibrium mentioned earlier. We show 
that the change in the momentum distribution as a function of 
modulating frequency to, Sn^u), carries information about 
the interacting spectral function both in the Mott and the su- 
perfluid side of the phase diagram. Using a non-perturbative 
strong coupling expansion|5l, we compute <5rik(w) and pro- 
vide a precise algorithm to extract spectral information from 
the measurement. We demonstrate that in the Mott phase, the 
variation of dn^u) in the Brillouin zone with u> tracks the 
dispersion of the gapped particle-hole excitations while in the 
superfluid phase, Sn^ shows a gapless linear dispersion at low 
uj and a gapped amplitude mode at high u. This leads to a 
qualitative distinction between the two phases. The height of 
the peaks of Sn^ in the Brillouin zone also contain important 
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information about the spectral weights of the different modes. 
Note that our proposed method contains detailed momentum 
resolved information which is not accessible by looking at the 
total energy transferred by lattice modulation suggested ear- 
lier 1231 l24l . Recently the study of these spectral features 
through Bragg spectroscopy (which measures the structure 
factor) has been suggested and the corresponding response 
functions have been calculated 11241 . However, compared to 
Bragg spectroscopy, our method is simpler to implement as it 
does not require additional lasers and the data can be collected 
in a single shot. Additionally, our method can yield momen- 
tum resolved spectral information, while Bragg spectroscopy 
integrates over the single particle momenta. Thus our pro- 
posal constitutes a significant improvement over the existing 
experimental proposals for both accurate detection of the SI 
transition and for obtaining the momentum-resolved spectral 
information for bosonic excitations in ultracold atom systems. 

We consider an ultracold atomic implementation of a re- 
pulsive Bose Hubbard model on a d-dimensional hypercubic 
optical lattice described by the Hamiltonian 



(a) 



(i) 



where b] creates a boson on the site i and n, 



b\bi is the 

number operator on site i, and we shall focus on d = 2, 3 in 
the present work. Here J, U and /i are respecively the tunnel- 
ing parameter, the onsite repulsion and the chemical potential 
of the system. At T — this system undergoes a SI transition 
as a function of J/U (or density or equivalently p/U) with 
the Mott insulating phase occuring at small J/U and at a com- 
mensurate integer density on each site. On top of this system, 
we consider a small perturbation to the optical lattice potential 
modulated at a frequency uj. Since the tunneling is exponen- 
tially dependent on the height of the optical lattice while other 
parameters have much weaker dependence, the primary effect 
of this perturbation is to modulate the tunneling parameter. 
The perturbing Hamiltonian can then be written as 



H 1 = -AJcos(wt) 2^6]6j, 

(ij) 



(2) 



where t denotes time and A is a small perturbation param- 
eter. For a deep lattice, J ~ E^exp{— \/Vq/Eii) where 
Vo is the lattice height and Er is the recoil energy. Hence 
^ ~ \/Vo/ ErSVo/Vq, where 5Vq is the amplitude of mod- 
ulation of the optical lattice height. This gives an order of 
magnitude for A. More quantitatively accurate estimates can 
be easily obtained by solving the band structure problem. 

We begin with the Mott phase of the model. Within linear 
response (linear in A), the change in n k is 



5n^{t) = rii(k, uS) cosujt + ri2(k, uj) sinwi, 



(3) 



where the in phase and out of phase response of the system, 
m and ri2, are given by the real and imaginary part of the 
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FIG. 1: The out of phase response of the momentum distribution in 
the n = 1 Mott phase of 2D Bose Hubbard model (U = 30 J and 
fj, = 20J) for (a)w = 20J and (b) uj = 35J. (c) and (d): Integrated 
response from 3D Bose Hubbard model (U = 40 J and p = 10.5 J) 
for (c) uj — 20 J and (b) uj = 45 J. The perturbation parameter 
A = 0.1 for all the plots. (Color online) 



response function n(k, uj + i0 + ) with 

n(k, iuu) = -p e k ^ G(k, ioj n )G{k, ioj + ioj n ). (4) 



Here /3 = l/T, T being the temperature of the system, 
ej. = —2 J X)i=i d cos ^» ' s tne b are band dispersion of the 
system and G(k, u>) is the single-particle Green function of 
the system. The response function in eq. |4] neglects vertex 
corrections although self-energy corrections can be included 
in the Green functions used to evaluate it. Note that for a non- 
interacting system n 2 is identically zero as the perturbation 
commutes with the non-interacting Hamiltonian. 

The single particle Green function in the Mott phase within 
the random phase approximation (RPA) 1 5 1 is given by 



G(k,iw) = z k (iuj - E+)- 1 + (1 - Zk)(iw 



(5) 
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where the particle (hole) dispersions are given by E k 

-8(i + | ^e k ± V e k + 4ekf a; + U 2 ^j and the spectral 

residue is given by z k = (E£ + dp, + Ux)/(E£ — E^). 
Here x = n + 1/2 and dp, = p, — U(n — 1 /2) where n is the 
integer number of bosons on each site. Using the above form 
of the boson Green function, we then obtain 



ni(k,oj) 
n 2 (k : uj) 



2Ae k 



A k z k (l - z k )[n B (E k ) - n B (E+)} 



7rAe k .z k (l - z k )[n B (_E k ) - n B (E£) 
x[<S(w + A k )-*(w-A k )], 



(6) 
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where Ak = E£ — EZ is the momentum dependent gap to 
creating particle-hole excitations in the system with zero cen- 
ter of mass momentum and tib is the Bose distribution func- 
tion. It is clear from the form of Ak that it has the lowest (but 
finite) value Am at the zone center (k = [0,0]) and increases 
as we go out towards the edge of the Brillouin zone. The SI 
transition is then marked by the vanishing of Am 

We focus here on U2(k, ui), which shows sharp features in 
the response. For modulation frequencies below the Mott gap, 
n 2 (k,ui) is zero everywhere in the Brillouin zone. For fre- 
quencies above the Mott gap n 2 will have response at the con- 
tours corresponding to Ak = to with a weight proportional 
to £k(l — ^k). Looking at ni(k, ui), we see that for frequen- 
cies less than the Mott gap, the response is broad with a broad 
peak at the zone center. For frequencies above the Mott gap, 
the response is sharply peaked around the contours of Ak = u> 
and changes sign as this contour is crossed. This qualitative 
reasoning is verified in Fig. [TJa) and (b), where n 2 (k,uj) is 
plotted as a function of k x and k y for d = 2 and two differ- 
ent frequencies corresponding to (a) Am < u> < U and (b) 
bj > U. For frequencies above the Mott gap, we see a contour 
of excitations which moves from the zone center to the edge 
of the Brillouin zone as the frequency is increased. Thus, one 
can determine the momentum dependent gap by doing a fre- 
quency sweep and looking at the peaks in the response. We 
note that for the 3D Bose-Hubbard model, the measured mo- 
mentum distribution is integrated along one axis and hence an 
analogous plot will reflect J_ dk z n 2 (k, u>)/2tt. The column 
integration leads to a region of excitations with the perimeter 
defined by the contour of excitations in the k x — k y plane. This 
is plotted in Fig.[T](c) and (d), for two different values of u. 

We now turn our attention to the superfluid phase, which is 
characterized by the presence of a macroscopic condensate 
and finite expectation values of the anomalous propagators 
F(k,iw„) = {bjj>\). In this case, the response function 
IT(k, ioS) in Eq. Elhas an extra term coming from the anoma- 
lous propagators and is given by 

n(k,io; n ) = ~4ek ^ \G(k, iu n )G(k, iu + iui n ) 

ioj n 

+ F(k,ioj n )F(k,iuj-iuj n )]. (7) 

In the superfluid phase, the normal and the anomalous Green 
functions are given by ||5) 



G(k, iu! n 
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52k 



F(k,iw n ) = $ 



iui n — Eik 
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icj n — E 2 k 
h 
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g4k 
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(8) 



B k ± V B l ~ Ck with Bk = b k - a 2 j2 + 

i 

^ + zJ, and bk = S^ 2 - C/ 2 /4 - (e k 



where ^ (2)k 

$ 2 /2andC k = 6^-$ 2 (^ + C/) 2 . Here a k = 2<5^-e k -2$ 
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FIG. 2: (a) and (b) : The out of phase response of the momentum 
distribution in the superfluid phase for a 2D Bose Hubbard model 
(U = 5J and fj, = 2.3J) for (a) u> = 4J (b) uj = 7 J. In (a) the 
response is due to creation of two phase modes. In (b) the low k re- 
sponse is from simultaneous excitation of a phase and an amplitude 
mode. There is also a small response from excitation of two phase 
modes at large k. (c) and (d): Column integrated out of phase re- 
sponse for 3D Bose Hubbard model (U = 8 J and n = 2.2 J) for (c) 
ui = 4.2J (b) ui — 11. 2,1 . The rings seen in the 2D response are 
now filled up due to column integration. The perturbation parameter 
A = 0.1 for all plots. (Color online) 



2$)(/i + U), with z — 2d for a hypercubic lattice in d di- 
mensions. The quasiparticle weights are given by gi( 2 )k 



[Eik ±(j* + U)][E 2 k T Ok^ik + bk}/(2Eik[E 2 n 



E 2k 



]). 



fl , 3(4)k is obtained from gu2)k. by taking Eik <-> E 2 k, and 
/k = (Efk - (P + U) 2 )/(E 2 U - £ 2 2 k ). The pole E 2k is the 
gapless phase fluctuation mode at small k (near the zone cen- 
ter). The pole Eik represents the amplitude mode at low k 
and is gapped. At large k the amplitude and phase modes are 
mixed; however, we still refer to these modes as amplitude 
and phase mode for easy identification. The phase mode car- 
ries almost all the spectral weight at low k, while the spectral 
weight resides mainly in the gapped mode at larger k values. 

Using these forms for the propagators we can now evaluate 
the momentum distribution response function. For simplic- 
ity, we restrict ourselves to T = and u> > (note that the 
in-phase and out of phase responses are respectively symmet- 
ric and antisymmetric in uS). Then, the change in momentum 
distribution is given by 



ni(k, oS) — 27rAek 



E 



Vk 



lk 



E 



2k 



uj 2 - (E 



-ak 



2Eik 
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E 2 k) 2 
2E 2 k 



n 2 (k,w) 



7rAe k [vk5(ui - E lk - E 2k ) 
+ctk5(uj - 2E lk ) + (3kS(u> - 
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where f k = -pikS4k - .92kff3k + $ 2 /k(l - /k)/(2-Bik^2k), 
«k = -5ikff2k + $ 2 /k/(4S 1 2 k ), and /3 k = -g 3 k54k + $ 2 (l- 
/k) 2 /(4£ 2 k ). 

Eq.[lO]shows that the out of phase response in the SF phase 
is peaked at energies corresponding to the creation of (a) two 
phase modes (b) an amplitude and a phase mode and (c) two 
amplitude modes. The in-phase response also shows broader 
peaks at these contours. At low frequencies (less than the gap 
in -Bik), the experiments should see a single contour of excita- 
tions corresponding to two phase mode excitations, as shown 
in Fig |2|a) for d = 2. As the frequency crosses the gap of 
the amplitude mode, two contours should be visible, a large 
contour corresponding to two phase modes and a small k con- 
tour corresponding to an amplitude and a phase mode, as in 
Fig |2|b). As the frequency crosses twice the gap, the two 
phase modes disappear as the frequency is above their band- 
width. The large k response is due to one amplitude and one 
phase mode. There is a low k response from two amplitude 
modes, but it is heavily suppressed due to spectral weights. 
For d = 3, a similar pattern is observed. However, the con- 
tours are replaced by regions with boundaries given by the 
location of the corresponding excitations in the k x — k y plane, 
as shown in Fig.|2jc) and (d). 

Thus we see that the change in momentum distribution due 
to optical lattice modulation can be used to obtain the exci- 
tation spectrum in both the Mott and the superfiuid phase. In 
the Mott phase, the particle-hole spectrum is gapped and there 
is no low frequency response. For both d = 2 and 3 and for 
modulation frequencies above the Mott gap, the response is 
peaked at the contour in the Brillouin zone whose dispersion 
matches the modulation frequency. Thus the peak response 
occurs on contours of increasing size as the frequency is in- 
creased. On the other hand, in the superfiuid phase, the exis- 
tence of gapless phase modes leads to low frequency response 
on contours whose size increase with frequency. However, 
on crossing the gap for the amplitude mode, there is a large 
response near the zone center and two distinct contours can 
be seen in this regime. At higher frequencies, there is a third 
contour of excitations corresponding to two amplitude modes, 
where the response is suppressed due to low spectral weight. 

In conclusion, we have derived the response of the momen- 
tum distribution of a Bose Hubbard model to optical lattice 
modulation within a strong coupling expansion, both in the 
Mott insulator and the superfiuid phase. In real experiments, 
the sharp poles in the single particle Green functions will be 
smeared due to several factors leading to a broadening of the 
response shown. The harmonic trap will smear the distribu- 
tion on a momentum scale of the order of the inverse oscil- 
lator length scale. For a wide trap, this is usually negligi- 
ble compared to other sources of smearing. Finite tempera- 
ture will lead to a smearing on energy scales ~ T. For the 
Mott insulating phase, this smearing would be negligible as 
long as the temperature is much smaller than the Mott gap 
(T <C If) 1 25 26 1, which is easily achievable experimentally. 
For the superfiuid phase, the low energy excitations would be 
smeared and it would be impossible to follow the very long 



wavelength low energy response (below energy ~ T). How- 
ever, the strong momentum dependence of the spectral weight 
of the phase modes would still lead to relatively well defined 
peaks in the Brillouin zone. On the other hand, deep in the su- 
perfiuid phase, the amplitude modes are gapped and the emer- 
gence of the small contour around the gap frequency should 
be easy to observe. Finally, including vertex corrections on 
top of the random phase approximation would broaden the 
poles, especially for the amplitude modes. The broadening 
represents processes where the amplitude mode can decay into 
two or more phase modes. However, for lattice systems, the 
phase modes have a finite bandwidth and a large gap for the 
amplitude mode restricts the phase space for such processes. 
Further, the two phonon decay involves phonon modes at rel- 
atively large energies which carry very little spectral weight. 
Thus the amplitude modes would be sharply defined near the 
zone center and hence the emergence of the small response 
contour around the gap should be easily observable. 
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